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Abstract 

The apparent universality of Bloch's T 3 / 2 -law for the temperature 
dependence of the spontaneous magnetisation, and of generalizations 
thereof, is considered. It is argued that in the derivation one should not 
only consider the exchange interaction between the spins, but also the 
other interactions between them, leading to elliptical spin precession 
and deviations from the parabolic dispersion of magnons. Also interac- 
tion effects are important to explain the apparent universality of gener- 
alized Bloch law exponents es, defined by M S (T) — M s (0)— const. -T tB , 
valid in a wide temperature range T\ < T < T2, and for dimension- 
alities d = 1, 2, and 3. The above-mentioned temperature range, the 
'Bloch range', lies above the quantum range, where magnetic long-range 
order (e.g. in d = 2 dimensions) is nontrivially enforced by the addi- 
tional interactions, but below the thermal critical region, where univer- 
sal 'anomalous scaling dimensions' apply. In contrast, for the Bloch 
temperature region, the universality is only apparent, i.e. a crossover- 
phenomenon, and simple scaling considerations with 'normal dimen- 
sions' apply. However, due to interactions, the Bloch exponent e B 
depends not only on the dimensionality d of the system, but also on 
the spin quantum number s (mod(l/2)) of the system, i.e. for given d 
the Bloch exponent es is different for half- integer s and for integer s. 
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05.50 Fh Phase Transitions: General Studies ; 
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1 Introduction: 

This paper grew out of discussions following a recent presentation 
of certain experimental results in our institute, which apparently, and 
unexpectedly, showed the T 3 / 2 Bloch law for M S (T) for a nanostruc- 
tured planar system, [1]; and although I originally thought that the 
considerations presented below were too simple for publication, some 
of the participants of those discussions suggested that I should write 
them down. So here I do so, just hoping to broaden and intensify 
thereby more personal discussions at the expense of being necessarily 
rather informal. 

Bloch's famous T 3 / 2 -law, and the natural generalization thereof, is 
written as follows: 

M S (T) = M s (0) - const. • T tB , (1) 

where according to the original derivation of F. Bloch for a d — 3- 
dimensional Heisenberg model, [2], the Bloch exponent €b would be 
3/2, while a recent analysis of various experimental data in different 
dimensions by U. Kobler and coworkers, see Ref. [3] and references 
therein, led to a set of values to be presented below, which are partially 
much different from the value 3/2, but appear to be quite universal. 
In particular, the number e# = 3/2 does not only follow from the 
original theoretical analysis of F. Bloch and the later refinement of this 
analysis by F.J. Dyson, [4], but (unexpectedly) it was also observed by 
U. Gradmann et al. in d = 2-dimensional ultrathin ferromagnetic films, 
see e.g. Ref. [5], and now it is apparently also observed in the planar 
nanostructures of Ref. [1]. Moreover, the experimental analysis of U. 
Kobler et al, Ref. [3], led to the above-mentioned set of apparently 
universal numbers for the Bloch exponent e#, where t B depends not 
only on the dimensionality d of the ferromagnetic or antiferromagnetic 
system considered, but also assumes different values for integer rsp. 
half-integer spin quantum number s: namely e# = 5/2 rsp. =3/2 rsp. 

2 for d = 1 rsp. = 2 rsp. 3 and half-integer s; whereas e# = 3 rsp. 
= 2 rsp. 9/2 for d = 1 rsp. = 2 rsp. 3 and integer s. (The e#-values 
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for d — 1 rsp. = 2 are also observed for anisotropically disordered or 
amorphous d = 2- rsp. = 3-dimensional systems.) 

Of course, as a theorist, one should always be skeptical against ex- 
perimental data which may be biased by theoretical claims. Therefore, 
in the present context it is stressed right from the beginning that the 
following presentation is not based on these data, but only motivated 
by it. 

In particular, and intentionally, I always write 'apparently universal' 
instead of simply 'universal', since the Bloch behaviour of Eqn. (1) is 
generally neither observed in the region of extremely low temperatures 
near T = 0, where quantum behaviour dominates, nor is it observed in 
the genuine thermal critical region, where one would obtain scaling with 
so-called anomalous dimensions, i.e. the 'critical exponents'. Instead, 
the above-mentioned behaviour should be considered as a crossover- 
phenomenon, observed in a (wide) intermediate temperature region 
Ti < T < T 2 , see below. In particular, the scaling dimensions involved 
correspond to simple powers of the basic thermal energy ksT rsp. of 
(/cbT) 1 / 2 , where ks is Boltzmann's constant. But in agreement with 
the fact that magnon-magnon interactions (although not yet thermally 
critical) already plays a role in the Bloch region (particularly at d — 2), 
the Bloch exponent, according to the considerations presented below, 
depends on the spin quantum number s in the above-mentioned way. 

The schedule of the present paper is as follows: In the next chapter 
the original Bloch derivation for a d = 3-dimensional Heisenberg system 
is considered, and, for later heuristical purposes, a simplified 'poor 
man's' variant thereof. Then I consider planar systems (d = 2). Next, 
d — 1 and d = 3 follow, all three cases for s = 1/2. Finally I consider 
s — 1, thereby collecting the whole set of 'apparent Bloch universality 
classes' of Kobler. At the end, in an additional chapter, the conclusions 
are presented. 

2 Bloch's law for d = 3 

According to this law (which is derived in any textbook on solid-state 
magnetism and is one of the most prominent results of theoretical 
physics) the temperature dependence of the magnetisation M S (T) of 
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a three-dimensional ferromagnet is given by 

M S (T) = M s (0) - const, x f d 3 k l — , (2) 

J exp (3e(k) — 1 



BZ 



because each excited magnon reduces the magnetic moment of a ferro- 
magnetic sample by a Bohr magneton. Here T is the Kelvin tempera- 
ture, f3 = (ksT)~ 1 with the Boltzmann constant k B , and the fraction 
{exp(/9e(/c)) — represents the thermal expectation value of the num- 
ber of magnons with excitation energy e(k), where the wave- vector k 
has its usual meaning. The integration is over the Brillouin zone BZ 
of the crystal. Since for k ■ a <C 1 (where a is the lattice constant) the 
excitation energy has a parabolic dispersion, e(k) oc k 2 , the integral in 
Eq. (2) converges for d = 3, whereas for d — 2 it would be 'infrared 
divergent' in agreement with the famous Mermin- Wagner theorem, [6]. 

For the excitation energy e(k) of a magnon in a Heisenberg fer- 
romagnet (and also in itinerant ferromagnets when the spin-orbit in- 
teraction, and also the Stoner excitations, are neglected against the 
collective magnon-like spin excitations) one simply has the following 
dispersion (we assume cubic symmetry and extend the validity of the 
formula over the whole integration region): e{k) = D • k 2 , where D 
is the so-called spin-wave stiffness. Therefore, by the substitution 
x := (3D ■ k 2 and the replacement d 3 k = 47ik 2 dk one gets the famous 
result M S (T) = M s (0) - const'. • (^) 3/2 , where (up to exponentially- 
small terms) the constants const, and const.' are related as follows: 

const.' = const, x 2tt f — ^. The analysis has been refined in the 

J cxp X— 1 J 

monumental work of F.J. Dyson, see Ref. [4], but this is not important 
below, except at the end. 

Rather, for Anderson's "poor man" , [7], instead of the usual deriva- 
tion, one can give the following simpler argument, which later-on will 
also serve for heuristical purposes: exp f3Dk 2 — 1 is approximated for 
long enough wavelengths and/or high enough T by the 'quasi-classical 
thermal-energy approximation' j3D-k 2 , so that one simply gets M S (T) = 
M s (0) - const. (^) x PhsR(T), where the factor k B T represents the 
thermal equipartition energy, while PhsR(T) means a typical phase- 
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't/o' 2 

space radius in k-space 1 , replacing the integral / ^prd/c x gk(T), 

where gk(T) is a temperature-dependent dimensionless weight-function. 
Here 'ir/a' represents a (very large) wavenumber-cutoff corresponding 
to the upper edge of the Brillouin zone, which is replaced by a sphere 
as in typical renormalization group arguments. But it would be wrong, 
if at this place one would perform directly the integration without any 
idea of the weight-function (after having made the above-mentioned 
'equipartition approximation' leading to the 'quasi-classical thermal- 
energy prefactor' (^fp))- Instead, one gets the correct PhsR(T) by 
a simple 'Pippard argument', i.e. simply by equating the dominating 
energy- resp. temperature-ranges: e(/c dom ) — D • k\ ora = ksT, i.e. 
PhsR{T) = k dom .(T) (- Ak(T)) = (¥)V2. 



3 Planar systems 

What would be different in a planar system magnetised in-plane ? 
One would again expect a quasi-classical 'equipartition approximation' 
exp (3e(k) — 1 — > (3e(k), but otherwise one would have a lot of differences 
and, generally, quite subtle behaviour: 

• First, instead of the simple parabolic dispersion e(k) = D ■ k 2 , 
representing excitations with circular spin precession, one would 

have a more general 'square root' formula, e(k) = \J e a {k) ■ e&(A;), 
representing an elliptical precession with the main axes a and b, 
respectively. 

For example, if one is dealing with a film of infinite extension in 
the x- and y-directions, with almost neglegible thickness th, then 
(if the film is magnetised in-plane, e.g. in the x-direction, see 
below) spin-wave deviations in the z-direction are strongly dis- 
favoured due to the demagnetising field H® M = —4tiM z , where 
M and M z are the z-components of the demagnetising field and 

— * 

the magnetisation, respectively. Thus, one would have e^k) = 
Cb + D ■ k 2 (= Cb in the long- wavelength limit, where Cb = 
Hb ■ 4:nM s is the effective anisotropy energy corresponding to the 

1 i.e. PhsR(T) is some kind of reciprocal thermal de Broglie wavelength 



5 



demagnetising field (/i£ is Bohr's elementary magnetic moment, 
the 'magneton'.)). In contrast, spin deviations in the y-direction 
would not be disfavoured, i.e. e a (k) = D ■ k 2 , as before (but see 
the following text, which should make things more precise): 

Actually, the magnetostatic interactions giving rise to the above- 
mentioned demagnetising field — 4irM z , also produce in-plane field 
components, which are proportional to the (small) thickness th of 
the ultrathin film. (Here, 'ultrathin' means that th is of the order- 
of-magnitude of some lattice constants; thus, even a monolayer 
should be approximated in this way by a continuous film of finite 
thickness th = a, and one should not simply set th = 0, instead.) 

• Second, the in-plane magnetostatic fields are not constant, but 
in contrast to the exchange interactions and to the usual single- 
ion anisotropies, they are long-ranged and induce a finite in-plane 
value of the magnetisation, although for the present planar sym- 
metry the Mermin- Wagner theorem, see Ref. [6], in the presence of 
only short-range interactions would rigorously predict vanishing 
long-range order. Actually, the above-mentioned finite magneti- 
sation is not in contradiction with this theorem nor with a recent 
extension, [8] , since in these papers the magnetostatic interaction 
is not included. 

• Third, generically, by the magnetostatic fields and/or other inter- 
actions, one has an effective out- of -plane anisotropy, favouring 

1. an out-of-plane state of the magnetisation (this case will be 
considered later), or 

2. an in-plane state, a case, on which we concentrate first: 

In the above-mentioned situation, for vanishing thickness, i.e. in 
the formal limit th — > 0, this in-plane state would be described by 
a model of short-range XY-symmetry, which has a phase transi- 
tion to a low-temperature phase with diverging susceptibility, but 
without long-range order for any finite T, the 'Kosterlitz-Thouless 
transition' (KT transition), see Ref. [9]. But for finite thickness, 
according to an equation to be presented below, there is always 
a /c-dependent long-range effective in-plane uniaxial anisotropy 
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favouring e.g. a spin-alignment along the (arbitrary) ±x-direction 
and disfavouring spin fluctuations with finite in-plane k along the 
±2/-axis, i.e. perpendicular to the magnetisation. (E.g. accord- 
ing to Hoffmann's well-known 'ripple theory', [10], the so-called 
'transverse ripple', with k = ±(0, k, 0), is disfavoured with respect 
to the longitudinal one, where k = ±(/c,0,0), see also Ref. [11].) 

— * 

• In fact, if the magnitude k of the (in-plane) wavevector k times the 
thickness th of the film is sufficiently low, due to long-range mag- 
netostatic interactions, for in-plane magnetised thin films with 
purely dipolar and exchange interactions, the magnon energy e(k) 
is given by the following formula, see Ref. [12], Eqs. (11)— (14): 

(3) 

where the dots denote neglected terms oc t\ or higher powers 
of th', i.e. to avoid z-dependencies it is assumed in the following 
that the thickness th of the film, although > 0, is nonetheless 
much smaller than the magnetostatic exchange length l e defined 
below. Furthermore, \i B is the Bohr 'magneton', /i B ■ jj- = D 
the spin-wave stiffness, M s the saturation magnetisation in c.g.s. 
units, and H x a small external field in the x-direction. (As men- 
tioned above, for th — this Zeeman field would be necessary to 
break the XF-symmetry which for short-range interactions would 
imply vanishing long-range order. Actually, the same magnetisa- 
tion divergence giving rise to the term oc k^/k 2 in Eqn. (3) (note 
the non-analytic behaviour on the direction!) would also give a 
contribution oc k y k x /k 2 modifying H x : this is usually neglected. 
However, in principle the neglection is not allowed for small H x ; 
e.g. the additional terms give rise to a 'blocking phenomenon' of 
the magnetisation ripple, see Ref. [10]. But here we neglect this 
addition as usual, since we do not consider that phenomenon.) 

• At the same time, due to the fact that by the strong elliptic- 
ity of the spin-precession the in-plane amplitudes are extraordi- 
narily large compared with the out-of plane amplitudes, which 
reduces the expectation value of the magnetisation much more 
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than in case of circular precession, the integrand on the r.h.s. of 
Eqn. (2) is multiplied by an additional ' ellipticity factor ' f(k) = 
{e a (k)+e b {k))j2 ^ w j 1 j ( j :i j s a i wa y S > i- j_ e _ now one would get instead 

V e a (k)-e b (k) 

of Eqn. (2): 

'tt/o' 

M S (T) = M S (0)- const, x / d 2 A; — * - • f(k) . (4) 

^ exp(pe(/c)) — 1 



o 



(Below, this still seemingly harmless and natural-looking equation, which 
is obtained by a standard Bogoliubov-Valatin transformation, see e.g. 
p. 25 of [13], from the bilinearized (i.e. approximate !) spin Hamilto- 
nian of the system, will be modified by spin interactions. This leads to 
a further temperature dependent phenomenological 'interaction factor' 
Fk(T), which is important in the Bloch region of d = 2-dimensional 
systems.) 

If one would again consider the formal limit t h — > of the Eqs. 
(3) and (4), then the integral on the r.h.s. of Eq. (4) would diverge, 

— * 

although the fraction in front of f(k) is now no longer dangerous. But 
the 'ellipticity factor' f(k) itself diverges: f(k) — > k~ l for k ■ th — > 0, 
and this ensures the logarithmic divergency of the integral in Eqn. (4). 
In fact, this divergency reflects the short-range XY-symmetry in the 
limit th = and is in agreement with the Mermin- Wagner theorem, 
see [6], since now all interactions have become short-ranged. But since 

for finite, but very small thickness, one has e(k) oc (k ■ th) 1 / 2 Ip-] 1 ^ 2 , 
a careful integration shows not only that the integral (4) converges 

'tt/o' 2tt 

(the dangerous part is oc f f ^4^4%- — ixit, which is ac- 

V to 1 k J =Q (k-t h )(iiTM a ) sm 2 (}!. fc + |^fc 2 ' 

tually very subtle since both t h and also A/M s must be different from 
0), but according to [14] one also gets a T 3 / 2 -behaviour at very low 
temperatures. However this is not relevant, since according to [14] 
this T 3 / 2 -behaviour is only observed in the lowest temperature-range, 
where the dipolar coupling is all-important, e.g. < T < T M • [y), 
where Tm = f^B^M s / ks is a characteristic temperature for the dipolar 
'shape anisotropy', while l e = \J 2 ^m 2 ^ s ^ ne aheady- mentioned 'mag- 
netic exchange length', [15]. 
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In Permalloy, with AtiM s = 10 4 Gauss and A = 1.3 10~ 6 erg/cm 
(remember Cb = [±b 4:TcM s and D := fi B (2A/M S )), the temperature 
T M would correspond to 0.67 K. As mentioned above, the temperature 
region for the Block behaviour is, instead, rather high: T\ < T < T2. 
Estimates for the two crossover temperatures Ti and T 2 are: 

• Ti = y/2 ■ he 4:TcM/kB] this is slightly larger than T M , but still 
very low compared to the temperature corresponding to the max- 
imal spin-wave frequency, T max : = D • (^) 2 /k B (= 845 K, with a 
lattice constant of a = n 1CT 8 cm, say). 

• But it is known that Block's T z l 2 -law applies to very high temper- 
atures, e.g. up to T 2 ps T max /3 ; without significant corrections, in 
planar systems of the above-mentioned kind, [5, 3]. And it should 
be stressed again at this place that this wide temperature region 
(the 'Bloch region') is still below the thermal critical region, [16]. 

In contrast, the relation e(k) oc (k ■ t h ) 1 / 2 Ipr] 1 / 2 is only valid for 
k ~ k c , where k c := th/l 2 . is a small 'crossover value': For larger 

wavenumbers the exchange interaction starts to dominate, i.e. f(k) — > 
1, and 

'tt/o' 

M,(T) » MAT ) - const x / 2^ ^' ^ - ; , (5) 

where k B T := e(k c ). By numerical evaluation (see below) this would 
inevitably lead to some kind of approximate 'Doring behaviour', [17], 
i.e. to M S {T) = M s (f) - const.'" x T ■ In f, where f can be considered 
as a temperature unit entering e.g. by quantum effects enforcing long- 
range order (In fact it turns out below that T is identical with Ti). 
As far as the author knows, the 'Doring behaviour' has never been 
observed, although, in the temperature range considered, it was also 
suggested by Maleev, [14], in a paper preceding that of Bruno, [18]. 

• Instead, to get the observed T 3 / 2 -behaviour (see [5, 3]) also in the 
'Doring temperature range', i.e. also forT > T\, one should note 
that, as soon as the wave-number k is large enough, there is no 
longer any in-plane uniaxial anisotropy felt by the magnons. 
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This means that it is natural to assume, as in the d = 2-dimensional 
XY model, that there exists a strong enhancement of the destructive 
character of these magnetic excitations, 'destructive' with respect to 
the magnetic long range order, leading to an additional temperature- 
dependent factor Fk(T) which represents the interaction of the mag- 
netic quasiparticle excitations. In the planar XY model, this additional 
'destructive factor' acts on all A;-scales (T-scales) until k — (T = 0); 
therefore, the magnetic long range order is completely destroyed in the 
XY-model. Physically, this happens through the dynamic formation 
of magnetic vortex excitations: In the low-temperature phase of the 
KT transition these vortices are strongly bound in pairs of opposite 
topological charge (opposite chirality), and the typical distance of the 
vortex centers of such 'vortex pairs', (1)t, remains finite, if the critical 
temperature is attained: At this temperature there is vortex unbinding 
at (l) Tc = ly. 

• As a consequence, spinwaves are not influenced by the vortex for- 
mation, if their wavelength A = 2ir/k is much larger than ly, 
whereas spatial spin-spin correlations with distances much smaller 
than ly are wiped out by the vortices. 

Similarly, we assume that in the Bloch region there is for k > 
k c a 'reduction factor' Fk(T) representing the enhanced reduction of 
magnetic order by the interactions between magnetic excitations of 
wavenumber k. In the present case, this reduction factor is not as 
fatal as in the d = 2-XY model, since wavenumbers smaller than k c 
are not involved. But one has again some kind of Pippard heuris- 
tics: Namely, along the lines of the original arguments of Kosterlitz, 
[19], we assume that F k (T) oc (k) T ■ W, where (...)t denotes a ther- 
mal average, so that the reduction increases with increasing T, i.e. 
increasing (&)t- In fact, according to the (by now already known) 
Pippard argument, (k)x can be replaced by (ksT / D) 1 ^ 2 . Thus again 
M S (T) = M s (f) - const, x T x T 1 / 2 , where on the r.h.s. the factor T 
represents the 'equipartition energy', while (for d = 2) the factor T 1 / 2 
arises from the 'reduction factor': This factor represents the strong 
reduction of long-range order by spin interaction effects, which are 
strongly 'relevant' in d = 2. (In d — 3 dimensions the 'reduction by 
interaction' is apparently not so important in the Bloch region. There, 
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instead, the magnetic moment of whole columns of sites with identical 
planar coordinates (x,y) is reduced by the additional fluctuations cor- 
responding to dynamic noncollinearities of the spins in the z-direction. 
This effect for d = 3 is already described with noninteracting magnons 
by our previous factor PhsR(T). Thus, d — 2 is much more dan- 
gerous to interactions than d = 3, in agreement with the fact that 
the 3-dimensional XF-model is 'harmless' compared with the d = 2- 
dimensional XY- problem.) 

There is another argument concerning Fk(T) oc T 1 / 2 : Since the typ- 
ical thermal energy of an entity representing one degree of freedom is 
ksT, and since the interaction should be dominated by two-body scat- 
tering, one guesses that in systems with strong magnon-magnon inter- 
actions the quasiparticle numbers (n^T — {exp(f3e(k) — should be 
multiplied in the present context by (/c^T) 1 / 2 : This concerns the first 
particle of an interacting pair. A similar interaction factor (ksT) 1 ^ 2 
should come from the second particle; i.e. together the two interacting 
particles might be considered as an entity involving the genuine ther- 
mal energy ksT = [(/c^T) 1 / 2 ] 2 . This is actually again some kind of 
Pippard argument; it will be repeated later, when the transition from 
s = 1/2 to s = 1 is considered in section 8. 

4 Quantum effects and crossover-flow for d=2 

It is important to understand what is going on: Therefore in the follow- 
ing we sketch a flow-line scenario for an in-plane magnetised thin film. 
Generally, we consider in the following a cartesian coordinate system 
where on the abscissa axis ('x-axis') we plot the effective out-of-plane 
anisotropy energy C& (e.g. C b = /is47rM s ), which favours an in-plane 
magnetisation and is at the same time proportional to an effective 
uniaxial anisotropy energy (e.g. e(k) = hb^kM s \{k ■ th/2) sin 2 y^] 1 / 2 ) 
favouring an (arbitrary) ±a>direction as in-plane magnetisation direc- 
tion. 

On the ordinate ('y'-axis) of our plot the reduced Kelvin tempera- 
ture ksT is presented. The origin of the cartesian coordinate system 
is at (C g = 0,T = 0). 

On the third axis (Vaxis) the magnetisation M s (C b , T) is presented, 
or any other thermal expectation value depending on C5 and T. 
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Along the x-axis for T = 0, one must consider quantum fluctuations 
("zero-point fluctuations"), which by the ellipticity of the spin preces- 
sion reduce the magnetisation from the equivalent of the spin quantum 
number s to an effective value, which is smaller: 

{k\Sf\k) = s - (k\Sf\k)=s-\ l2 (k)\ 2 . 

Here \k) := £ -7=- ex P( iA; * rr ) \GS), where \GS) is the (non- 



/2sN 

n 



trivial !) ground state of the system. The vectors / count the N sites 
ff of the lattice considered, the S-* etc. are the spin operators. 

By the Bogoliubov-Valatin transformation, see e.g. page 25 in Ref. 
[13], the ground-state \GS) is not identical with the 'trivial state' |0), 
where all spins are aligned in the preferred direction (i.e. |0)=|toj = s) 
V/, i.e. for all sites of the lattice). Instead, \GS) is obtained from 
|0) by a straightforward calculation, which hides the 'ellipticy' by a 
mathematical formalism described in detail in the above-mentioned 
'text-book reference'. In this way one derives the following identity 
involving the 'ellipticity factor' f k = M fc )+j*(*0)/ 2 f Eqn. (4) and the 

'zero-point spin reduction' |72(fc)| 



2. 



\i2(k)\ 2 = \-(h-i)= (6) 



-1) 



1^ {2A/M s )k 2 + 2ttM s + ttM s ■ {k ■ t h ) ■ sin 2 n 

2 ^((2A/M s )k 2 + 4ttM s ) • ((2A/M s )k 2 + 2vrM s • (k ■ t h ) ■ sin 2 

In Fig. 1 the zero-point spin reduction |72(^)| 2 is plotted as a func- 
tion of k for sin 2 ip^ = 1 with Permalloy material parameters (i.e. 
4ttM s = 10 4 Gauss, A = 1.3 x 1CT 6 erg/cm) and a thickness of th — 4 
nm, as in the magnetic nanostructures studied by Hollinger et al, [15]. 

As can be seen from Fig. 1, |^y 2 ((0, k, 0))| 2 is less than 10 % for 
k > 0.1 (nm)" 1 ; for these A;- values the linearized Bogoliubov-Valatin 
theory makes sense. Whereas for ^-values below that value, corre- 
sponding to almost macroscopic wavelengths, there is a strong increase 
of the ellipticity so that for these small wavenumbers the linearization 
would be no longer meaningful. 



But for th = one sees that for k-l e 1, where l e := = ^ 

is the above-mentioned exchange-length (ss 5.7 nm for Permalloy), a 
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divergency develops for k — > 0: |72(^)| 2 ~ \ • f% — l/(4/c • / e ) ^> 
1. This divergency has already been mentioned above, and in the 
present context for T = it means that for small enough values of 
k ■ l e spatially periodic spin fluctuations of this small wavenumber k 
cannot be treated any longer in the Bogoliubov-Valatin theory (i.e. 
by linearization), since e.g. the effective spin-quantum number s eS = 
s — \l2(k)\ 2 becomes negative for these k. This is true although the 
integral / d 2 k {^(k)] 2 is finite (whereas, as already mentioned above, 
for th = and finite temperatures, long-range order would not at all 
exist, i.e. J d 2 k \tl2&g±l = In contrast, for t h > 0, |7 2 (A;)| 2 would 

diverge in a similar way, namely as 0.25/ \J (k ■ th/2) sin 2 tp^ + k 2 ■ I 2 , 
where the first term depends essentially only on the thickness and not 
on the other parameters, but long-range order at finite T exists in this 
case nonetheless for all values of So the convergence of the integral 
for M S (T) and the behaviour of the zero-point reduction \^2{k)\ 2 for 
small k are not simply related. 

But as a consequence of the preceding paragraphs, due to the mag- 
netic dipole interaction there is automatically a k-dependent effective 
in-plane uniaxial anisotropy favouring the axis parallel to the mag- 
netisation. This k-dependent effective anisotropy can be roughly ob- 
tained by replacing the effective field 2nM s ■ k ■ th • sin 2 ipk by its angu- 
lar average 2ttM s ■ k • th • 0.5. Here the typical long-wavelength scale, 
e.g. the typical wavelength in the observation of magnetisation ripple 
phenomena, see Ref. [10], is given by the characteristic wavenumber 
k — l/l e . So we associate to the magnetostatic phenomena in thin 
films a k-dependent effective in-plane uniaxial anisotropy of the or- 
der of C u {k) := hb • nM s ■ k ■ th, i.e. we write C u {k) = C u • k, with 
C u := hb ■ nM s ■ t h - Here k < l/l e . 

In Fig. 2, again for the already mentioned Permalloy parameters, 
we present results for the 'naive magnetisation decay' M S (T), i.e. cal- 
culated according from the formula (4) for d = 2 with the 'ellipticity 
factor' f k , but without the 'interaction factor' F k (T). By comparison 
with similar results for d — 3, where the T 3 / 2 behaviour is quite evi- 
dent, we get for d = 2 the already mentioned 'Doring behaviour', i.e. 
the result looks like (M S (T) - M S (T )) oc T-ln £, where T is a temper- 
ature representing an effective 'cutoff' at low k (see the figure caption 
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of Fig. 2). 



Since the experimental observations deviate apparently from the Doring 
behaviour in a significant way, the additional 'interaction factor', men- 
tioned above, seems really necessary. 

On the horizontal coordinate axis of our plot one has thus an 'or- 
dered line segment' ranging from the above-mentioned value C u , cal- 
culated e.g. for k = \/l e or any other relevant characteristic recipro- 
cal length of the system, up to infinity, whereas the dangerous 'in- 
frared ellipticity fluctuations' corresponding to smaller wavenumbers 
are 'gapped away' at very low temperatures, since the amplitudes be- 



From the value on the x-axis corresponding to the above-mentioned 
standard case (e.g. Permalloy) on the 'ordered line segment' a flow 
line emerges, which is well-defined if along the third axis (V-axis) 
a quantity like M s (Cb,T) is plotted, and which flows to the thermal 
critical point. 

This line begins at T = 0, i.e. in the quantum regime, with a positive 
vertical slope, but soon it may turn to the left, since with increasing 
T the typical wavenumbers increase, the ellipticity decreases, and the 
exchange interactions become important. So the flow line in this region 
may resemble (details are not needed in the following !) a flat elliptical 
segment running almost parallel to the x-axis in the direction of the y- 
axis, until it extrapolates, almost horizontally, to the above-mentioned 
small crossover temperature Ti (= \^2Tm, see below). This should 
signal the crossover from the linear rsp. square root behaviour e(k) = 



Cb ■ Dk 2 + C u ■ k of the dispersion for k-l e < 1 to parabolic behaviour, 



e(k) ^ D-k 2 for k-l e > 1. Therefore one roughly expects k B T x ^ V2C b , 
see e.g. Eqn. (3) for th = 0: It is, in fact, visuable from Fig. 3, that 
there is a rather sharp crossover of this kind at k c := Titj-Jl 2 = 0.1 
(nm)^ 1 . At shorter wavelengths (larger k), the exchange interaction 
dominates and the interaction energy is ~ ksT\ or larger. 

It is not essential (and not at all universal) how this 'crossover point' 
is reached, i.e. from which ground state at T = (as we have seen 
the ground states will be nontrivial, i.e. not given by the Bogoliubov- 
Valatin transformation at length scales k • l e <S 1). In contrast, in 
the Bloch region Ti < T < T 2 , where the dispersion of the magnons, 
whether interacting or not, is dominated exclusively by the exchange 
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interaction, the behaviour is essentially universal, e.g. the values of the 
anisotropies, including the effects of the magnetostatic interactions, 
become irrelevant. 

Analogously, at the second axis ('y'-axis) one has the following sce- 
nario: long-range magnetic order up to the thermal critical temper- 
ature (Curie temperature) T c , which is of the order of the maximal 
magnon energy, T max := D ■ (^) 2 /k B , or typically even larger, [16]. 
From this Curie fixed point, i.e. in the thermal regime, our (backward) 
flow line emerges starting at first horizontally, but soon turning down- 
wards, until it flows back from T 2 : = T max /3 towards T\ essentially 
according to Bloch's T 3//2 -law, thereby reaching the above-mentioned 
crossover point, where ksT\ ~ V2Cb, as already noted. From there 
the (backward) flow returns to the starting-point (corresponding e.g. 
to ultrathin Permalloy films with finite th <C h) on the 'ordered line 
segment' of the where the linearized (i.e. Bogoliubov-Valatin) 

transformation should be valid until the very small effective in-plane 
uniaxial anisotropy of the order of C u = C b ■ ^: But this last part of 
the flow is non-universal, i.e. C u can also be replaced by a quantity not 
depending on Cb- 

Now, what counts in the Bloch behaviour between T\ and T 2 = 
l^max; ar & the thermal energy factor oc T and the effective interac- 
tion radius (in k-space) oc T 1 / 2 corresponding to only the segment of 
the flow line between T\ andT 2 (in this region the quasi- classical 'poor 
man's scaling' considerations apply). 

5 Two-dimensional systems with positive out-of- 
plane anisotropy 

For two-dimensional systems with a positive uniaxial perpendicular 
anisotropy, i.e. favouring (in contrast to the previous case) an out-of- 
plane orientation, there is no ellipticity of the spin precession around 
the z-direction, and the temperature dependence of the magnetisa- 
tion follows in a naive formulation analogously to Eqn. (4), namely 

M S (T) = M s (0) - const. • vr x /° gg* with e(k) = C 6 out + D ■ k\ 

k=0 

where C out := — C& — 4ttM s > 0, whereas in-plane magnetostatic fields 
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are now neglegible. 

The above-mentioned integral can be evaluated analytically, yield- 
ing a variant of the 'Doring behaviour', namely 

M S (T) = M s (0) - const. • « x f^V ( -£^ML + ln ^ ~ , 
sV ; sK J \ D J \ k B T e x ™ - 1) ' 

(7) 

with e max := /3 • (CT + D ■ (f ) 2 ) and e min := (3 • C 6 out - 

Again, if in accordance with [3] a T 3 / 2 -behaviour would be observed, 
an additional T 1,/2 -factor would be necessary. 

In the particular case of a 'reorientation transition', corresponding 
to a situation, where the effective uniaxial anisotropy would be zero, 
C£ ut = 0, for short-range interactions the original Mermin- Wagner the- 
orem would become applicable and the magnetic order would vanish: 
Instead, again by the (long-range) dipolar interactions, one observes a 
nontrivial domain structure. 

At this place it should be be mentioned that the 'reorientation tran- 
sition', and in particular spinwave excitations (including the dipolar 
interactions and surface anisotropies) in systems with a finite number 
of non-equivalent magnetic planes with different tilt angles of the mag- 
netization direction in these planes, have been carefully studied in a 
series of papers by Erickson and Mills, [20, 21, 22]. 



6 d=l: T 5 / 2 -behaviour ? 



Here we consider a homogeneous ferromagnetic wire stretching from 
x = — oo to x — +oo along the with constant circular cross sec- 

tion, magnetised longitudinally, i.e. in the x-direction. Now the dipolar 
(effective) anisotropy acts with equal strength, Cb = Hb • 2nM s , along 
the y-axis and the z-axis, so that, by applying the 'naive approach' of 
Eqn. (4) again, we get (almost as before) M S (T) = M s (0) — const, x 

7r/a 

J — d ^ , with ek = Hb ■ 2irM s + D-k 2 . For temperatures allowing 

q GXp jj ' €.yfS J _L 

the approximation exp((3e(k)) — 1 = (3e(k) this can be directly evalu- 
ated and yields M S (T) = M s (0) — const, x T x n/2y/fj,B ■ inA, whereas 
for general temperatures one gets the numerical behaviour presented 
in Fig. 4, i.e. also a behaviour oc T. 
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That this behaviour cannot be valid to arbitrary low temperatures 
follows from the Maxwellian relation dS/dH = dM/dT and Nernst's 
'Third Fundamental Theorem' of thermodynamics, according to which 
dS/dH -> for T -> 0. 

In fact, experimentally one has obtained instead, see [3], a pro- 
nounced T 5,/2 -dependence, in agrrement with the apparently universal 
classification of U. Kobler for half-integer s and one-dimensional sys- 
tems, [3], and astonishingly also with the experiments by U. Gradmann 
and coworkers, [23], on Fe monolayers on W(110) (covered by Ag). 
Also the films with the different T 3 / 2 -behaviour studied in [5] are pla- 
nar Fe films on W(110), but thicker (from 3 to 21 monolayers), and 
whereas (as already mentioned above) in [5] one has the T 3//2 behaviour 
discussed in the two preceding sections, for the films of [23] one has in- 
stead the T 5//2 behaviour, which is still more astonishing. Not only for 
this system, but also for d = 1 this is unexpected and suggests that here 
the 'interaction damage' factor F k (T) is even stronger, namely oc T 3 / 2 
instead of oc T 1 / 2 . Presently, this is not at all understood, or only very 
qualitatively: One knows that in one dimension the critical thermal 
fluctuations destroy long range order at finite T even for Ising models, 
whereas in d — 2 the fluctuations are substantial for XF-models, but 
not for Ising models. So it is not quite surprising that already for the 
Bloch region the interaction factor should be much stronger for d — 1 
than for d = 2. 

Topological excitations may play an essential role also in this case, 
namely the so-called Bloch-point excitations, which are studied in a 
recent paper of Thiaville et. al, [24]. These are very 'floppy'; they 
consist topologically of opposite spin directions along a certain axis, 
e.g. S oc ±z for z > rsp. z < 0, but in the xy-plane there can be 
a vortex. In this case, the 'damage factor' F k (T) corresponding to 
vortex interactions exists for all three orthogonal planes, the xy-plane, 
the yz-plane and the zx-plane. 

In any case, the T 5 / 2 -behaviour in the Bloch region does not con- 
tradict the fact that ultimately, i.e. after a further crossover to the 
genuine thermal critical region, there can exist not even Ising critical 
order in the wire, such that finally T c — > 0. For Gradmann's mono- 
layer films, [23], instead, after the pseudo-one-dimensional (non-Ising) 
behaviour in the Bloch region, the system should return to the genuine 
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two-dimensional critical (Ising-) behaviour near T c . 

So this is quite complicated and subtle, as usual if one is dealing 
with dipolar interactions. 

7 T 2 crossover-scaling for itinerant magnets in d=3 

The arguments of the following section look especially simple (although 
this may be deceptive): We consider a crystalline metallic ferromagnet 
in d=3 dimensions. The system has separate energy bands {E^(k)} 
rsp. {Ei(k)}, and at T = the 'Stoner excitations', i.e. electron-hole 
excitations with spin flip, where e.g. an electron with energy 'immedi- 
ately below the Fermi energy Ep\ with initial-state wave- vector k and 
spin |, is moved to a state 'immediately above Ep\ with a different 
final-state wave- vector k and spin [. Now at finite T, the expression 
'immediately above or below Ep y means: 'within an interval of width 
a* ksT around Ep. Every such excitation reduces the magnetic mo- 
ment of the sample by two Bohr magnetons; i.e. in the above-mentioned 
'heuristical formula' one has a 'phase-space factor' oc (/cb^) 2 , where one 
of the two powers of k B T stands for the effective phase-space of the ini- 
tial states, while the other one represents the final states. There is no 
'energy factor', since all this happens in the above-mentioned interval 
immediately near Ep. 

Similar arguments also apply to itinerant antiferromagnets, in ac- 
cordance with the observations of U. Kobler, [3] , who always observed 
the same Bloch exponent for the magnetisation of ferromagnets rsp. the 
sublattice magnetisation of antiferromagnets (As a consequence the 
Bloch scaling is probably not directly related to quantum-criticality, 
since the quantum critical behaviour of ferromagnetic and antiferro- 
magnetic systems is different). In contrast, for semiconducting sys- 
tems as EuO and EuS the magnetism should better be described 
by a (non-itinerant) d = 3 Heisenberg model with s = 7/2, and so 
one expects the T 3 / 2 -behaviour for d = 3. But an 'interaction fac- 
tor' of T 1 / 2 would transform this into the observed T 2 behaviour, in 
accordance with [3]. Yet such a behaviour would definitely go be- 
yond the Heisenberg model, since within that model there is the rig- 
orous analysis of Dyson, Ref . [4] , for arbitry fixed spin quantum num- 
ber s, see Ref. [4], predicting a more gradual decrease: M S (T) = 



18 



M s (0) - ci • T 3 / 2 - c 2 • T 5 / 2 - c 3 • T 7 / 2 — c 4 • T 4 — C(T 9 / 2 ). Note that 
the first term with integer exponent is oc T 4 , and not oc T 3 ; the integer 
term results mainly from the decrease of the magnon stiffness, which 
according to rigorous theories 'a la Dyson' should be oc T 5 / 2 , whereas 
more phenomenological theories, and according to Kobler, Ref . [3] , also 
the experiments, suggest a decrease following an effective behaviour as 
D(T) = D(0) — const. • T 3 / 2 . 

Because of the large value of s = 7/2 for Eu 2+ , the 'non-Dysonian' 
behaviour, and the 'extra' interaction factor oc T 1 / 2 , look not com- 
pletely unreasonable. And one should note again that this behaviour 
goes definitely beyond the (underlying) Heisenberg model, which al- 
ways assumes fixed s, in contrast to the Pippard arguments from the 
previous section and those discussed below in connection with the prob- 
lem of the difference between integer and half-integer s. And one should 
also note, concerning the case of metallic magnetism, that the Stoner 
mechanism (which is actually a two-particle mechanism) looks com- 
pletely different from anything formulated in the Heisenberg model. In 
contrast, the collective magnon excitations, which are also present in 
itinerant magnets and which have the same quadratic dispersion as be- 
fore, seem to give a neglegible contribution to the behaviour of M S (T) in 
the temperature range considered, although at very low temperatures a 
T 3 / 2 -contribution would always be larger than a T 2 -one. So this means 
again that one should remain above a certain crossover-temperature, 
in this time that one to the dominance of the T 3 / 2 behaviour for d = 3. 

8 Integer s 

The arguments in the previous sections have all considered half-integer 
values of s, e.g. s = 1/2 and s = 7/2. For integer values, e.g. s = 1, 
one must consider the additional interaction energy, ~ ksT, which is 
necessary to break e.g. a s = 1 triplet state of two s = 1/2 spins 
into two separate s — 1/2 entities. This happens at sufficiently high 
temperatures at the upper range of the Bloch region, i.e. below the 
genuine critical temperature region. 

So for a single entity the square root of that energy comes into 
play as an additional energy factor so that Kobler's effective Bloch 
exponents of = 3/2 (rsp. = 5/2) for spin systems with half-integer 
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s and (effective) Bloch-dimensionality d = 2 (rsp. d — 1) should change 
for integer s to e# = 4/2 (rsp. e# = 6/2), apparently as observed, 
see [3]. For d = 3 things are more complicated: according to the 
tentative arguments of the previous section, the Bloch exponent t B = 2 
for d — 3 and half-integer s is given by a product of three terms, 
which can be written symbolically as (k B T)} nitial statc x (k B T)\ nal statc x 
(k B T)°, where the third term represents the constant energy-factor. 
However for integer s, due to the above-mentioned interaction effects, 
the first two factors should each acquire an additional power of k B T, 
while the third factor should acquire an additional 'square-root power' 
(k B T) 1 / 2 . In this way, starting with the T 2 -behaviour for half-integer 
s one would expect a T 9/,2 -behaviour for integer s, as observed, [3]. 
Whether this is more than heuristics, has to be shown, similarly to the 
previous sections, but it should be noted that the pronounced crossover 
from a T 2 low-temperature behaviour to a T 9 / 2 behaviour at higher 
temperatures, which U. Kobler observes for NiO at T = 100 K rsp. for 
Fe at 628 K, seems to be in favour of the present scenario 2 . 

9 Conclusions 

As a consequence of the preceding sections, concerning the temperature 
range considered, it should be noted that the arguments of the present 
paper, supported by quantitative calculations especially in sections 2 
and 3, do neither concentrate on the magnetostatic interactions (al- 
though these are very important, as seen in those chapters) nor stress 
the immediate vicinity of the quantum region at T = nor that of the 
thermal critical region near T c , but rather concentrate on the (wide) 
'Bloch temperature range' between the (rather low) crossover temper- 
ature T\ (= \/2 Cb/k B ) and a (much higher) crossover temperature 
of the order of magnitude of T 2 = \ D ■ (n/a) 2 jk B . 

So quantum behaviour is significant, but at the same time, in the 
Bloch region, it only plays a secondary role, apparently in contrast 
to interactions, which appear to be particularly important for d = 2, 
where they are known to destroy the long-range order for the Kosterlitz- 

2 At several places in the observations of U. Kobler, see [3], also a crossover from a T 3 ^ 2 - 
behaviour 'at low temperatures', e.g. for T < 200 K, to a T 2 -behaviour at a considerable 
range of higher temperatures is observed. 
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Thouless phase transition. In higher dimensions also the main implicit 
assumption of the Heisenberg model, that the spin quantum-number s 
is invariably constant, becomes questionable through the interactions 
already in the Bloch region. 

All statements in the previous sections, except from section 8, apply 
to half-integer s in Kobler's classification, [3]. Already at this place, i.e. 
for half-integer s, the role of the interactions seems crucial, particularly 
for d = 2 and d—1. Additional complications for integer s (outside the 
region of overwhelming importance of the thermal fluctuations) have 
been discussed in section 8. 

A more thorough analysis of the questions left open, which mainly 
concern the role of interactions, is shifted to future work. 
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Figures and Figure Captions 




Fig.l: The zero-point spin reduction |72(A:)| 2 (see the text, especially 

— * 

Eqn. (3)) is presented over the wavevector k = (0, k, 0) , where k 
is measured in units of (nm)" 1 , for an ultrathin Permalloy film of 
thickness th = 4 nm, in a double-logarithmic plot with material pa- 
rameters 4nM s = 10 4 Gauss, where M s is the cgs-magnetisation, and 
A = 1.3 x 10~ 6 erg/cm the spin- wave stiffness. 

The characteristic magnetostatic exchange length l e = (2A/(47rM s )) 1 / 2 
is 5.7 nm, such that the characteristic value k c := th/{2l1) is ~ 0.06 
(nm)" 1 ; |72(&)| 2 vanishes at k — ir/l e i.e. at ps 0.5 (nm)" 1 . For 
k > 0.1 (nm)" 1 it is less than 10 %. For these A;-values the linearized 
Bogoliubov-Valatin theory makes sense. Whereas for /c-values below 
that value, corresponding to almost macroscopic wavelengths, there is 
a strong increase of the ellipticity so that for these small wavenumbers 
the linearization would be no longer meaningful. 
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Fig. 2: Here for dimensionalities d = 2 (left panel) and d = 3 (right 
panel), respectively, the integral AM S (T) := M s (0) — M S (T) oc 

/■ d d A; 
fc £ expjT- 1 y/[2 • (1 - cos k) + C b ] ■ [2 • (1 - cos fc) + fc • C u • 2 -§-}} - 1 

2-(l-cosfc) + ^ + A;-C , M -{f 

yj[2 ■ (1 - cosfc) + C 6 ] • [2 • (1 - cos A;) + k ■ C u ■ 2 -§-} 

is presented over the reduced Kelvin temperature T. 

This plot represents schematically (in the left panel) the temper- 
ature dependence of the magnetisation of an in-plane magnetised ul- 
trathin film under the influence of a typical exchange interaction e = 
2 • (1 — cos A;), which behaves for small k as e = k 2 , while the reduced 
dipolar anisotropy oc 4itM s is represented by C& = 0.01, and, finally, 
a /c-dependent effective in-plane uniaxial anisotropy corresponding to 
the last term in Eqn. (3), e(Jfe) oc ^ ... + (4vrM s ) 2 • (k ■ t h /2) ■ {k 2 /k 2 ), is 

represented by the terms oc k ■ C u , where C u = 0.001 is chosen. The 
results of the left panel correspond roughly to the 'Doring behaviour' 
70 • T • In j£r=, which necessitates for d = 2 an additional interaction 

0.05 ' 

factor F fc (T) oc T 1 / 2 to obtain a T 3 / 2 -behaA r iour, whereas the 'd = 3- 
results' of the right panel (from which the difference can be seen) would 
yield just this behaviour. 
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Fig. 3: In a double-logarithmic plot the energy e(k) (up to a constant 
factor) is presented over k for k — (0, k, 0), for material parameters of 
permalloy as in the preceding figure, where k is measured in (nm)^ 1 . 
The thicknesses of the ultrathin films considered, are 4 nm, 2 nm, and 
0.6 nm, respectively, corresponding to the highest, medium rsp. the 
lowest curve. 

Note the sharp crossover around k = k c = itth/ll = 0.1 nmr x . For 
shorter wavelengths (larger k) the different thicknesses can no longer 
be distinguished, since the exchange interaction dominates. 
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0.2 0.4 



Fig. 4: For C b = 0.01, as in Fig.2, the integral M s (0) - M S (T) oc 

fc=7T 



/ 

fc=0 



dA; 



expjT- 1 ! 2 • (1 - cos A;) + C b \ - 1} 



is presented over the reduced Kelvin temperature T. The results should 
be valid for d — 1. 
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